In the Claims: 

1. (Original) A method for determining a geomagnetic transmission function for the 
transmission of a population of particles, said particles having one or more rigidities, to a 
satellite in a known earth orbit, under known geomagnetic conditions, comprising the steps: 

(A) dividing said orbit into a plurality of sequential steps comprising at least a first step 
and a last step, each sequential pair of said steps being connected at a point, each of said points 
corresponding to a position of said satellite along said orbit, and specifying a value for Kp and a 
value for Dst to specify said geomagnetic conditions; 

(B) for a plurality of said points, selecting a plurality of arrival directions for the arrival 
of said particles at said satellites, wherein each of said arrival directions represents either an 
allowed or a forbidden trajectory of one of said particles; and 

(C) for each of said arrival directions at each of said plurality of said points, determining 
whether said arrival direction represents an allowed or a forbidden trajectory by tracing a path of 
said particle to arrive at said point with said arrival direction, until said particle path intersects a 
boundary of earth's magnetosphere, thereby indicating an allowed trajectory, or until said 
particle path intersects either earth's atmosphere or earth's surface, thereby indicating a 
forbidden trajectory, wherein said tracing of said path is performed by integrating in the time 
domain the Lorentz equation 




wherein F is the force vector acting on said particle, m is mass, y is 1/(1 -v 2 /c 2 )' 1 ' 2 , v is said 
particle's velocity vector, t is the travel time of said particle, Q is the charge of said particle, and 
B is the magnetic field vector acting on said particle, wherein B is given by 



B = BioRF(r s t')+ BTSYo(Kp,r,tO 



wherein B IG rf is the International Geomagnetic Reference Field promulgated by the 
International Association of Geomagnetism and Aeronomy, and B TS yo is the modified 
Tsyganenko field given by the sum of B^™ + B Ymg ro + + B^ 0 ^ + By^^ + 

B^ C) + BxJV + B Ysm (C) + B^ + Bxsm^ + B Y J™ + BzJ™, wherein coordinates in the 
solar magnetospheric system are denoted sm, and wherein coordinates in the solar magnetic 
coordinate system are denoted mg; and wherein 

B x m = Q t xz r ; 

B Y ro = Q t yz f ; 
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wherein z, = z-Zs(x,y,^), 
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wherein S TRC = Jp 2 + (a TRC + , 

#T,RC = V 2 ?" 1 " ^T,RC 

D r = d 0 + # 2 + r AW+ rA(*)» 

wherein C, = -98.72 when Kp = 0, 0 + , Q = -35.64 when Kp = 1",1,1 + , Ci = 
-77.45 when Kp = 2',2,2 + , C,= -70. 12 when Kp = 3~,3,3 + , Q = - 162.5 when Kp 
= 4 _ ,4,4 + , Q = - 128.4 when Kp * 5", 

wherein C 2 = -10014 when Kp = 0, 0 + , C 2 = -12800 when Kp = C 2 = 

- 14588 when Kp = 2~,2,2 + , C 2 = - 16125 when Kp = 3\3,3 + , C 2 = - 15806 when 
Kp = 4",4,4 + , C 2 = - 16184 when Kp * 5", 

wherein a T = 13.55 when Kp = 0, 0 + , a T = 13.81 when Kp = r,l,l + , a T = 15.08 
when Kp = 2",2,2 + , a-f= 15.63 when Kp = 3~,3,3 + , a T = 16. 1 1 when Kp = 4",4,4 + , 
a T = 15.85 whenKp ^ 5", 

wherein D 0 = 2.08 when Kp = 0, 0 + , D 0 = 1.664 when Kp = r,l,l + , D 0 = 1.541 
when Kp = 2",2,2 + , D 0 = 0.9351 when Kp = 3",3,3 + , D 0 = 0.7677 when Kp = 
4",4,4 + , D 0 = 0.3325 when Kp ;> 5", 

wherein Rc = 9.084 when Kp = 0, 0 + , Rc = 9.238 when Kp = 1",1,1 + , R* = 9.609, 
when Kp = 2\2,2 + , R*= 8.573 when Kp = 3",3,3 + , = 10.06 when Kp = 4-,4,4 + , 
Rc= 10.47 whenKp* 5", 
wherein Ly = 10 R E , 

B Y (RC) = QRcyzr; 
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By (M) = e^tdoyz cos y + (C n y+C 12 y 3 +Ci 3 yz 2 ) sin y], and 

^(M) = e x/Ax [(Ci4 + Cisy 2 + CieZ 2 ) CQS y + ( CnZ+ Cl8Zy 2 + Ci9Z 3) sin ^ 

wherein C 6 through C 19 are given by: 
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and wherein Ly=10.0, D x =13.0, L RC =5.0, Lf=6.30, y^4.0, S=0.010, y^l.O, Rt=30.0, 

Xqc-4.0, L^SO.O, andD yc =20.0, 
and thereby determining whether said particle's trajectory intersects either the boundary of 
earth's magnetosphere, thereby indicating an allowed trajectory, or intersecting earth's surface or 
earth's atmosphere, thereby indicating a forbidden trajectory. 



2. (Original) The method of claim 1, wherein said plurality of arrival directions comprises 
one or more randomly or pseudorandomly selected arrival directions. 

3. (Original) The method of claim 1, further comprising the step: 

(D) for each particle rigidity, determine what fraction of particles of that rigidity will be 
transmitted to said satellite. 

4. (Original) The method of claim 1, wherein said plurality of points comprises a complete 
set of points for said orbit. 

5. (Original) The method of claim 1, wherein said plurality of points comprises a set of 
points for a portion of said orbit. 



6. (Currently amended) A method for determining, for a given particle environment outside 
of earth's magnetosphere, what portion of a population of particles having one or more rigidities 
making up said particle environment will be transmitted to a satellite in a known earth orbit, 
comprising the steps: 

(A) performing steps (A) through (C) of claim 1, thereby computing a geomagnetic 
transmission function for said population of particles; and 

(B) applying said geomagnetic transmission function to said population of particles, 
thereby determining what portion of that population of particles will be transmitted to said 
satellite. 

7. (Currently amended) A method for determining, for a given particle environment outside 
of earth's magnetosphere, what portion of a population of particles having one or more rigidities 
making up said particle environment will be transmitted to a satellite in earth orbit, comprising 
the steps: 

(A) prompting a user to specify an earth orbit; 

(B) determining whether said orbit is among a group of preselected orbits, each of said 
orbits in said preselected group of orbits having an associated predetermined geomagnetic 
transmission function for a range of particle rigidities, each of said predetermined geomagnetic 
transmission functions having been prepared in accordance with claim 1; and 

(G) for the case wherein said orbit is among said preselected group of orbits, applying 
said predetermined geomagnetic transmission function for said orbit to said particle 
environment outside earth's magnetosphere. 

8. (Original) The method of claim 7, wherein said preselected group of orbits comprises a 
quiet shuttle orbit having a 450 km altitude and a 28.5° inclination, a disturbed shuttle orbit 
having a 450 km altitude and a 28.5° inclination, a quiet space station orbit having a 450 km 
altitude and a 51.6° inclination, and a disturbed space station orbit having a 450 km altitude and 
a 51.6° inclination, wherein said quiet shuttle orbit has a geomagnetic transmission function 
given by FIG. 4, wherein said disturbed shuttle orbit has a geomagnetic transmission function 



given by FIG. 5, wherein said quiet space station orbit has a geomagnetic transmission function 
given by FIG. 6, and wherein said disturbed space station orbit has a geomagnetic transmission 
function given by FIG. 7. 

9. (Currently amended) The method of claim 7, further comprising the step: 

(B) for the case wherein said orbit is not among said preselected group of orbits, 
performing a step selected from the group consisting of (a) returning an error message to a user, 
and (b) computing a geomagnetic transmission function for said orbit, in accordance with the 
method of claim 1. 

10. (Original) A method for determining a flux of a species of solar ions having a specified 
atomic number between 3 and 92, and a specified kinetic energy, for a satellite in a near earth 
orbit, comprising the steps: 

(A) specifying an atomic number for solar ions for evaluation; 

(B) specifying a kinetic energy for said solar ions; 

(C) specifying a baseline model for said flux of solar ions, wherein said baseline model is 
selected from the group consisting of a worst day model, a worst week model, and a peak flux 
model; 

(D) in the case wherein said specified atomic number is greater than 20, selecting iron as 
an elemental spectrum model; 

(E) in the case wherein said specified atomic number is less than or equal to 20, selecting 
oxygen as an elemental spectrum model; 

(F) looking up a value for an elemental breakpoint, wherein said elemental breakpoint is 
a function of said elemental spectrum model and said baseline model, and wherein said 
elemental breakpoint is selected from the table: 



^""^^^^ Baseline model 
Elemental ^^^^^^ 
spectrum model ^^"^^^^ 


worst day or peak flux 


worst week 


iron 


24.23 MeV/nuc 


19.90 MeV/nuc 


oxygen 


15.94 MeV/nuc 


12.89 MeV/nuc 



(G) in the case wherein said baseline model is said worst week model, and said elemental 
spectrum model is iron, and said kinetic energy is greater than 127.93 MeV/nuc, calculating an 
unsealed flux, wherein said unsealed flux equals A™ x (En/ Me % uc ) wherein = 
3.16814xl0 6 (cm 2 sr MeV/nuc)' 1 , (En/ Me % uc ) is said kinetic energy, normalized to be 
dimensionless, and Gsp = 2.861; 

(H) in the case wherein the condition recited in step (G) is not satisfied, and wherein said 
kinetic energy is greater than said elemental breakpoint, calculating an unsealed flux wherein 
said unsealed flux equals A 3 xEN rJ1 , wherein A 3 is a function of said elemental spectrum model 
and said baseline model, and wherein said A 3 is selected from the table: 



Baseline model 
Elemental ^"^-^^^ 
spectrum model ^^^^^ 


worst day or peak flux 
in (cm 2 sr MeV/nuc)" 1 


worst week 

in (cm 2 sr MeV/nuc)" 1 


iron 


0.252948xl0 10 


0.249719xl0 9 


oxygen 


0.106702xl0 10 


0.667628xl0 9 



and wherein y si is a spectral index and is a function of said elemental spectrum model and said 
baseline model, and wherein said y si is selected from the table: 



^^^^^^ Baseline model 
Elemental 

spectrum model ^^^^^ 


worst day or peak flux 


worst week 


iron 


-4.52970 


-3.7610 


oxygen 


-4.14060 


-3.76850 



(I) in the case wherein the condition recited in step (G) is not satisfied, and wherein said 
kinetic energy is less than or equal to said elemental breakpoint, calculating an unsealed flux 
wherein said unsealed flux equals A 2 exp(-GxEn 1/4 )xEn 1/4 , wherein A 2 is a function of said 
elemental spectrum model and said baseline model, and wherein said A 2 is selected from the 
table: 



Baseline model 
Elemental ^^"""^^ 
spectrum model 


worst day or peak flux 
in (cm 2 sr MeV/nuc)" 1 


worst week 

in (cm 2 sr MeV/nuc)" 1 


iron 


1.8991xl0 8 


3.0372xl0 8 


oxygen 


4.9518xl0 8 


1.1307xl0 9 



and wherein G S j is a spectral index and is a function of said elemental spectrum model and said 
baseline model, and wherein said G S j is selected from the table: 



^^^^^^ Baseline model 
Elemental ^^^-^^^ 
spectrum model ^^^^^^ 


worst day or peak flux 


worst week 


iron 


5.70 


5.70 


oxygen 


5.70 


5.70 



; and 



(J) in the case where said atomic number is between 3 and 92, inclusive, calculating a 
solar ion flux for said specified atomic number and said specified kinetic energy by multiplying 
said unsealed flux by a scale factor ratio, said scale factor ratio being the ratio of a scale factor 
for an element having said selected atomic number over a scale factor for said spectrum model 
element, wherein said scale factors are selected from the table: 



Atomic 
Number 


ocale 
Factor 


Atomic 
Number 


ocale 
Factor 


Atomic 
Number 


ocaie 
Factor 


Atomic 
Number 


L>cale 
Factor 


Atomic 
Number 


Scale 
Factor 


Atomic 
Number 


Scale 
Factor 


5 


0 


22 


4.377 

vlft" 3 


39 


4.878 

vlO" 6 

X 1 V 


56 


4.878 

vlfl" 6 
X IvJ 


73 


2.195 

X 1U 


90 


4.878 


6 


4.704 
xlO' 1 


23 


4.088 


40 


1.22 
xlO' 5 


57 


4.878 
xlO* 7 


74 


2.439 
xlO" 7 


91 


0 


7 


1.205 
9 

xlO" 1 


24 


1.65 
xlO' 2 


41 


9.756 
xlO" 7 


58 


1.22 
xlO" 6 


75 


4.878 
xlO* 8 


92 


2.927 
xlO" 8 


8 


1 


25 


5.625 
xlO" 3 


42 


4.878 
xlO -6 


59 


1.951 
xlO" 7 


76 


7.317 
xlO" 7 


3 


0 


9 


4.560 

976 

xlO* 5 


26 


1 


43 


0 


60 


9.756 
xlO' 7 


77 


7.317 
xlO -7 


4 


o 


10 


2.131 
2 

xlO* 1 


27 


L303 
xlO' 2 


44 


2.195 
xlO -6 


61 


0 


78 


1.463 
xKT 6 






11 


1.744 

715 

xlO" 2 


28 


3.172 
xlO' 2 


45 


4.878 
xlO -7 


62 


2.439 
xlO" 7 


79 


2.439 
xlO* 7 






12 


2.062 
4 

xlO' 1 


29 


3.048 
xlO" 4 


46 


1.463 
xlO"* 


63 


9.756 
xlO" 8 


80 


2.439 
xlO" 7 






13 


1.826 

829 

xlO" 2 


30 


7.457 
xlO" 4 


47 


4.878 
xlO* 7 


64 


4.878 
xlO -7 


81 


2.195 
xlO* 7 






14 


3.593 
5 

xlO" 1 


31 


4.878 
xlO" 5 


48 


1.707 
xKT 6 


65 


7.317 
xlO" 8 


82 


2.439 
xKT 6 






15 


2.279 

675 

xlO" 4 


32 


1.22 
xlO -4 


49 


2.195 
xlO* 7 


66 


4.878 
xlO* 7 


83 


1.463 
xlO" 7 






16 


9.758 
xlO* 2 


33 


7.317 
xlO" 6 


50 


4.878 
xlO* 6 


67 


9.756 
xlO* 8 


84 


0 
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17 


1.680 

488 

xlO -4 


34 


7.317 
xlO" 5 


51 


3.415 
xlO* 7 


68 


2.439 
xlO -7 


85 


0 




18 


1.771 

545 

xlO' 3 


35 


9.756 
xlO^ 


52 


7.317 
xlO* 


69 


4.878 
xlO" 8 


86 


0 




19 


3.644 

715 

xlO* 4 


36 


4.878 
xlO" 5 


53 


1.463 
xlO -6 


70 


1.951 
xlO" 7 


87 


0 




20 


4.826 
xlO -2 


37 


7.317 
xlO -6 


54 


6.585 
xlO* 


71 


4.878 
xlO" 8 


88 


0 




21 


2.929 
xlO" 4 


38 


2.439 
xlO" 5 


55 


4.878 
xlO" 7 


72 


1.951 
xlO" 7 


89 


0 





wherein elements having atomic number of 20 or less are scaled to oxygen, and elements having 
higher atomic numbers are scaled to iron. 

11. (Original) A method for determining a flux of solar ions of species having a range of 
atomic numbers between 3 and 92, wherein each of said species includes particles having a 
range of specified kinetic energies, for a satellite in a near earth orbit, comprising the steps: 

(A) for each kinetic energy of each of said species, performing steps (A) through (K) of 
claim 10, and storing a flux value for each of said kinetic energies for each of said species; and 

(B) adding each of said flux values to obtain a total solar ion flux value. 

12. (Original) A method for modeling the effects of cosmic rays on microelectronics on a 
computer connected to the internet, comprising the steps: 

(A) receiving a login message from a remote user connected to the internet, and 
generating and transmitting back to said user with a script in response thereto an HTML main 
menu page, wherein said main menu page comprises prompts for said user to run one or more 
routines selected from the group consisting of: (1) a routine for calculating a geomagnetic 
transmission function for SEU-inducing particles, (2) a routine for calculating a flux of SEU- 
inducing particles in the near earth environment or in the environment shielded by earth's 
magnetosphere, (3) a routine for calculating a solid shielding transport function for SEU- 
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inducing particles, (4) a routine for calculating a proton induced single event upset rate, (5) a 
routine for calculating a linear energy transfer rate, and (6) a routine for calculating a heavy ion 
induced single event upset rate, said main menu page further comprises prompts for said user to 
select a user request file for each of said routines; 

(B) in response to a message from said user to run at least one of said routines, using a 
user request file specified by said user, generating and executing with a script a system command 
for said computer to run each of said specified routines; and 

(C) in response to the completion of all of said routines specified by said user, generating 
and transmitting to said user with a script an HTML page notifying said user of said completion. 

13. (Original) The method of claim 12, wherein said main menu page generated and 
transmitted to said user further comprises prompts for said user to create or edit a user request 
file for any one of said routines, in lieu of running said one or more of said routines, and further 
comprising the step: 

(D) in response to a message from said user to create or edit a user request file, 
generating and transmitting to said user with a script an HTML page with prompts for the entry 
of fields for said user request file; and 

(E) in response to said user transmitting information for said user request file, storing 
said information in a user request file with a script. 

14. (Original) The method of claim 12, wherein said main menu page comprises prompts for 
said user to run a routine for calculating a geomagnetic transmission function for SEU-inducing 
particles, and wherein said step of generating and executing with a script a system command for 
said computer to run each of said specified routines comprises at least generating and executing 
a system command for said computer to run said routine for calculating a geomagnetic 
transmission function for SEU-inducing particles, in accordance with claim 1. 
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15. (Original) The method of claim 12, wherein said main menu page comprises prompts for 
said user to run a routine for calculating a flux of SEU-inducing particles in the near earth 
environment or in the environment shielded by earth's magnetosphere, and wherein said step of 
generating and executing with a script a system command for said computer to run each of said 
specified routines comprises at least generating and executing a system command for said 
computer to run said routine for calculating said flux of SEU-inducing particles in the near earth 
environment or in the environment shielded by earth's magnetosphere, wherein said routine 
comprises calculating a flux of solar heavy ions in accordance with claim 10. 



